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ABSTRACT 



We use a simple perturbation theory argument and measurements of 
charmonium leptonic widths r{ipi^s ~^ e^e~) to estimate the ratio 
^. -Ro = |^r;cis(0)l /l^»/'is(0)l i^ ^^^ general context of non- relativistic po- 

r-| ! tential models. We obtain Ro = 1.4 ± 0.1. We then apply well known 

potential model formulas, which include lowest order QCD corrections, to 
k> ; find r(77c — *■ 77)/r(^is — >• e"'"e~) ~ 2.2 ± 0.2. The central value for 

j^ ' r(V'is ~^ e"*"e~)in the 1992 Particle Data Tables then leads to a (non rel- 

ativistic) prediction T{ric -^ 77) ~ 11.8 ±0.8 keV. This prediction is in good 
agreement with a recent measurement by the ARGUS collaboration, is con- 
sistent with a recent measurement by the L3 collaboration but is significantly 
higher than several earlier measurements and than previous theoretical esti- 
mates, which usually assume Ro = 1. The correction to i?o = 1 is estimated 
to be smaller but nonnegligible for the bb system. Using the current cen- 
tral measurement for r(Ti5 —>■ e^e^) we find T{rib —>■ 77) ^ 0.58 ± 0.03 
keV. A rough estimate of relativistic corrections reduces the expected two 
photon rates to about 8.8 keV and 0.52 keV for the rjc and rjf, mesons respec- 
tively. Such corrections however, are not expected to significantly affect our 
estimates of Ro. 



Potential models (non-relativistic [NR] as well as " relativized" versions) have been 
successfully used to describe many properties of quarkonium {cc and bb) states |^. Rel- 
ativistic effects and beyond lowest order QCD corrections are expected to be more im- 
portant in cc meson than in bb meson but it is hard to devise model independent tests 
(i.e. that do not depend strongly on the particular form of the potential being used) 
that may pinpoint some properties of charmonium spectroscopy or decays where a N-R 
description clearly fails. 

In the present paper we make a prediction for T^rjc ~^ 77) {^{Vb "^ 77)) which relies 
solely on i) a N-R description of the cc (bb) system, ii) experimental data of leptonic 
charmonium (bottomium) decays, iii) approximate validity of lowest order perturbation 
theory for the color-hyperfine splitting interaction and iv) approximate validity of lowest 
order QCD radiative corrections. At the end of the paper we make an estimate of 
relativistic corrections to our results. 

As a starting point we make use of a previously derived result |^ relying on a N-R 
description of quarkonium systems which includes lowest order QCD corrections: 



r«.is^e+e-) 3'- ■ 7r'(2m,,r|*^„,(0)|' 

where a^ should be evaluated at the charm scale. We will use as{mc) ~ 0.28 ± 0.02 [Q. 
For consistency of the N-R description of the system and since the hyperfine splitting 
will be included only perturbatively, one should set 2mc ~ M^ §. 
Thus, one obtains 

r(r7c^77) _ 1 ^^ l^r?c 



r(^i5^e+e-) • |vl/^,,(0)| 
and analogously, using as{mb) ^ 0.18 ± 0.01 



1-57 ,' '^-""„ , (2) 






0.37 ' "^ ^' 2 . (3) 



It is important to note that even if we were far more conservative with the uncertainty 
in the value of as(mc) and ^^(m;,), the resulting uncertainty in the numerical coefficients 
of eqns. (2) and (3) would only be of order of a few percent. 

A widely used approximation at this point is to set |\1'^^(0)| /|\l/^^g(0)| ~ 1 
(|vl>,jO)|V|vl/x,,(0)|' ^ 1) leading to T{r], ^ 77)/r(^i5 ^ e+e") ^ 1.6 
(r(?7b — > 77)/r(Ti5 — > e+e~) ~ 0.37). We will show below that for the cc case this 
commonly used assumption is off by more than 30% (and could be off by as much as 
50%). For the bb case the correction to this approximation is estimated to be smaller 
but still significant. 



In the context of N-R potential models, the interaction Hamiltonian responsible for 
the splitting between the l^Si states and the I^Sq state of a QQ meson is given by 
H = Hsi2 + Hso where 

^^12 = ~Z + y/^l-^Sf^ (^ = ^5i2 + ^5i2 • (4) 

/Ti and /r2 are respectively the static color- magnetic moment operator for the quark and 
the antiquark and r*is the relative position vector. Color indices have been omitted and 
electromagnetic interactions have been ignored. The Spin-Orbit part of the Hamiltonian, 
Hsoi is included for completeness. However, it does not contribute to lowest order to 
either the energy splitting (eqn. (5) below) or the shifts of the wavefunctions (eqns. 
(6-9) below) for S states. 

We now study the effects of the Hamiltonian H to lowest order in perturbation theory. 
The energy splitting between the l^Si and I^Sq states is given, to lowest order by 

AB<') = E%^ - B«„ =< *™ Ji/.J*™ . > - < *Sk|i/.J*Sk > 

= < *S°4.I^IJ*S°4 > - < *Skiff|.J*S?.. > . (5) 

where the last step follows after angular integration from the fact that the relevant states 
have £ = 0. Experimentally, AE^c = M^^^ - M^^^g = 118 ± 2 MeV @]. We estimate 
AE^ii = M-fis — M^^jg = 45 ± 15 MeV, using the measured value of AEcc and N-R 
potential model formulas which include lowest order QCD corrections. 

The lowest order correction to the radial and orbital ground state wavefunctions at 
a point r due to the Hamiltonian H is given (in common bra-ket notation) by 

,(0) 



A^^^^r^sAr) = A(^) < m,.s. >= E ' " \.(or",;o) ' ^^'^ ' (6) 



p(0) ,p(0) 



<f|vl/W ><M/W|/75 .J^S3k> 
i?i? - E^^ 



A«M/i.,„(f) = A«<rlM/i.,„>= E ' " ^(0) " ' (0) ° ■ (7) 

But the wavefunction at the origin is zero for all non-S states and therefore only 
S-states contribute to the sum for r = 0. Thus, the shifts in the wavefunctions at the 
origin are given as follows: 

A^'^issAO) = 1^ p(o) p(o) 

N>1 -C/15 - -C/TVS 






p{0) p{0) 



vi/(°)rnir < vi/W 1/75 ivi/W > 

(0) Z?(0)^iT,(0) 



^1 (i5i5^-i5;]v^^)^r3^5.(o) 

The last step is possible because the matrix element has support only at the origin. The 
previous step follows from the fact that only S-states are involved in the sum. Note 
that only one state per radial excitation contributes to the sum. The wavefunction in 
the denominator is a short hand notation to mean that multiplication on the left by 
< ^I'lS^jO > will cancel it. 



In the same way, 



,2 



Therefore, to lowest order in perturbation theory, 
1^135.(0)1' - l^i^sMf = [^i?5,(0)* (A«v[;,3,,(0)) + (A«v1/i3^^(0)*) ^^^i^) 



^%^{0)* (A«v[/,.,„(0)) + (A(^)v[/,.,„(0)*) M/f4(0) 

N>1 i^lS - ^NS 



2^^' z. 7^(0^ — ^^(or ' ^^°^ 



where Aii^'^-'^^ is given by eqn. (5). 
Thus, 



aWO)r _ i + aAgC)^ K(g)l' 



1*135,(0)1* "1 l*i»s,(o)|' 4°s - b!°s 



iV>l 1^135,(0)1 M^35, -M135, 

where the last step (which is consistent to lowest order) was taken to express the result 
in terms of experimentally measurable quantities. 

Notice that since we are working in the context of N-R potential models we can 
express the ratio of the wavefunctions at the origin for S vector states in terms of ratios 
of leptonic widths and masses [0] : 



I*iv35,(0)|' r(iV3Si ^ e+e-)M 



|v^i35,(0)|^ T{l^S,^e+e-)Ml,,^ 



N^S, (12) 



Our result is then 

We calculate this ratio using M^^^ — M^^^^g = 118 ±2 MeV as well as the information on 
masses and leptonic widths for the ^^i states ^(3097), ^(3685), ^'(3770), ^'(4040), ^(4160) 
and V^(4415) given in the 1992 Particle Data Tables @]. We obtain 

''^''^^^(°)''- 1.4 ±0.1 . (14) 



The error was estimated from the given experimental uncertainties 0] . We have 
ignored possible contributions from "continuum" states above DD threshold because 
these are outside the potential model Hilbert Space. We are assuming implicitly that 
the discrete spectrum constitutes a complete set of states of the relevant N-R Schrodinger 
equation and thus including (physical) continuum state contributions would amount to 
double-counting. 

We note that over half of the contribution of the sum over states in eqn. (13) to our 
result in eqn. (14) is due to the lowest excitation ?/'(3685) (~ 0.23) while the contribution 
of the highest observed radial excitation ?/'(4415) is quite small (~ 0.03). Therefore, we 
expect that the five radial excitations of ipis observed so far saturate the sum to a good 
approximation. 

Although we have no reliable estimate of the corrections to the result of eqn. (14) 
due to effects that are higher order in -^5^2) ^^^ magnitude of the lowest order correction 
(~ 40%) can be used as an indication that such corrections are not likely to change our 
result significantly (i.e. by more than 0(10%)). Such corrections should be viewed as a 
systematic uncertainty of our approximations and will be ignored hereafter. 

Combining our results in eqns. (2) and (14), we obtain the prediction 

(^c->77) _ 2.2 ±0.2 . (15) 



Using the value for the ipis leptonic width quoted in ref. 0], 

V{il) -^ e+e") = 5.36 ±lf^ keV, (16) 

we obtain an absolute estimate 

r(r7c -^ 77) = 11.8 ± 0.8 ± 0.6 keV, (17) 

where the first uncertainty stems from our main result eqn. (15), while the second 
reflects the experimental uncertainty in eqn. (16). 



We point out that if we had used the common approximation |\1/^^(0)| /|\&^(0)| = 1 
instead of our estimate (eqn. (14)), we would have obtained a rate T{ric -^ 77) ~ 8.4±0.5 
keV instead of our result in eqn. (17). This is in line with most previous theoretical 
predictions (around 7 keV)0 and also agrees within errors with the experimental av- 
erage (r(77c -^ 77) = 6.6 ±li keV) quoted in ref. g (See also refs. @, [0] and ||). 
A more recent measurement by the ARGUS Collaboration 0, on the other hand, is 
r(^c — ^ 77) = 12.6 ± 4.0 keV, which is centered closer to our prediction (eqn. (17)). 
Finally, a recently pubhshed result by L3 JlOl, ^{Vc -^ 77) = 8.0 ± 2.3 ± 2.4 keV, is cen- 



tered closer to the estimate based on the assumption Ro = I but is still consistent with 
our prediction (eqn. (17)). 

Although the effect is smaller for the bb mesons, we repeat the same procedure for 
completeness. Using our above mentioned estimate for the energy splitting AEi^i = 
M-fig — Af^j^jg = 45 ± 15 MeV and the experimental information about T states (masses 
and leptonic widths) 0, we obtain an estimate, analogous to the result in eqn. (14), of 

|2 

^ 1.16 ±0.06 . (18) 



l*...s(0)f 



Using our results, eqn. (3) and eqn. (18), we obtain 

n^t-^ii) 0.43^0.02 . (19) 

r(Ti5 ^ e+e ) 

The current experimental value for r(Ti5 -^ e^e") |^ then leads to the prediction 

^{Vb -^ 77) = 0.58 ± 0.03 ± 0.02 keV , (20) 

where the first error stems from eqn. (19) and the second one from the present experi- 
mental uncertainty in r(Ti5 -^ e+e^)[Q. 

To summarize, the main point of the present paper is to show that within the general 
context of N-R potential models one can (using available experimental information) 
make a reliable estimate for the ratio -Ro = |\E'ii5q(0)| /|\E'i35^(0)| for cc and bb systems. 
A commonly used approximation is -Ro = 1- Our estimates (eqns. (14) and (18)) imply 
that -Ro = 1 is a bad approximation for the cc system and requires moderate corrections 
for the bb system. 

The main measurable consequence of this study is that Tijjc — > 77) is expected to 
be significantly larger (by between 30%and50%) than predictions based on the -Ro = 1 
assumption. In fact, the N-R "potential model" predictions for the rates of all important 
decay modes of the ?7c (relative to r(?/' -^ e+e^)) are enhanced by the same amount with 
respect to the predictions based on the -Ro = 1 assumption. Thus, if -Ro = 1 is used for 
such predictions the error made would have to be compensated by using unphysically 
large values for the strong coupling as at tree level. 



From the theoretical viewpoint, precise measurement of T{r]c — > 77) favouring lower 
values (eg. around 7 keV, which agrees well with the central value in ref. [Q and also 
the recent measurement ref. [|10]) would put into question the validity of N-R potential 
model assumptions for the description of cc systems. The only way that we see to escape 
this conclusion would be to argue that the ratio (2;;^) iii eqn. (1) is significantly smaller 
than one. However, this seems to go against the weak-binding assumption that is needed 
for self-consistency of N-R potential models. 

On the other hand, if r(?7c — > 77) turns out to be close to our prediction (eqn. 
(17)) [See also the measured value in ref. 0] we could state that N-R potential model 
descriptions of charmonium have passed yet another test successfully. 

All the above results and discussions are based on a strictly N-R description of the 
cc and bb system. As stated at the beginning of the paper, relativistic corrections 
may be significant, especially for the cc system. We conclude our paper with a brief 
discussion of the likely effects of such corrections on our predictions. We think that 
our results for the ratios Ro (eqns. (14) and (18)) should be essentially unaffected by 
relativistic corrections, because the main input for those estimates consists of actual 
experimental data (leptonic widths and mass splittings) which of course include full 
relativistic corrections. The argument used for these estimates (eqns. (4) to (13)) relies 
more on lowest order perturbation theory than on strictly N-R dynamics. On the other 
hand, eqn. (1) and its analogue for bb systems does rely more directly on the N-R 
description. The main relativistic correction to that ratio comes from the fact that the 
contribution of the quark propagator between the two photons in the rjc {rjb) decays is 
sensitive to the momentum distribution of the quarks in the decaying meson [see, for 
TT| and [O ]. We roughly estimate these effects by using commonly 
=< /5^ > for the orbital and radial ground state cc and bb 
0.1. The /5-dependence in the integrands 
of eqns. (13) and (15) in ref. [ll[] is " pulled out " of the integrals by replacing /3^ by 
< /?^ >. This procedure results in a relativistic correction factor of 



example, ref. 

quoted values for ^| 

quarkonium: < jS"^ >^«i 0.25 and < /3^ >t 



l-</3^> 



In 



i+V<^ 



Vi-</3^ 



2 + Vl-</5'> 



to the RHS of eqn. (1) and its bb analogue. Using the above values for < (3^ > leads to 
a numerical factor of about j^ for the cc system and yu ^'^^ ^^^ ^^ system. Using our 
N-R central estimates given in equations (17) and (20) above, we obtain " relativistically 
corrected " central estimates of TuEiiVc —^ I1) ~ 8-8 keV and TuEbijIh -^ 77) ~ 0.52 
keV. 
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